##### Exerccios 4 - Anlises Exploratrias#####

#### 4.2 Cervejas ####

cervejas <- c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
cervejas <- factor(cervejas, levels = c("chope", "lata", "garrafa"))
library(lattice)
histogram(cervejas)
dotplot(cervejas, horizontal=F)


#### 4.3 Caixetais ####

cax = read.csv("caixeta.csv", header=TRUE, as.is=TRUE)
head(cax)
cax$dap <- (pi/4)*(cax$cap/10)
hist(cax$dap)

str(cax)
histogram(~h|local, data = cax)

#### No h diferenas representativas

##### 4.4 Eucaliptos ######

e.grandis <- read.table("egrandis.csv",as.is=F,header=T,sep=";")
e.grandis
boxplot(e.grandis$dap~e.grandis$regiao, ylab="dap")
boxplot(e.grandis$dap~e.grandis$rotacao, ylab="dap", xlab="rotacao")

str(e.grandis)
qqnorm(e.grandis$ht)
qqline(e.grandis$ht)


#### 4.5 Mais Caixetais  #####

cax <- read.csv("caixeta.csv", header=TRUE, as.is=TRUE)
cax$dap <- (pi/4)*(cax$cap/10)
unique(cax$local)
help(plot)
plot(cax$dap[c(which(cax$local=="chauas"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="chauas"),which(cax$especie=="Tabebuia cassinoides"))])
plot(cax$dap[c(which(cax$local=="jureia"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="jureia"),which(cax$especie=="Tabebuia cassinoides"))])
plot(cax$dap[c(which(cax$local=="retiro"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="retiro"),which(cax$especie=="Tabebuia cassinoides"))])

help("scatter.smooth")
scatter.smooth(cax$dap[c(which(cax$local=="chauas"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="chauas"),which(cax$especie=="Tabebuia cassinoides"))])
scatter.smooth(cax$dap[c(which(cax$local=="jureia"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="jureia"),which(cax$especie=="Tabebuia cassinoides"))])
scatter.smooth(cax$dap[c(which(cax$local=="retiro"),which(cax$especie=="Tabebuia cassinoides"))],cax$h[c(which(cax$local=="retiro"),which(cax$especie=="Tabebuia cassinoides"))])

help("xyplot")
xyplot(dap[(cax$especie=="Tabebuia cassinoides")]~h[which(cax$especie=="Tabebuia cassinoides")]|local,data = cax)
splom(cax[ , c("dap","h")],groups = cax$local[(cax$especie=="Tabebuia cassinoides")] )
histogram(dap[(cax$especie=="Tabebuia cassinoides")]~h[(cax$especie=="Tabebuia cassinoides")]|local,data = cax)
qqmath(dap[(cax$especie=="Tabebuia cassinoides")]~h[(cax$especie=="Tabebuia cassinoides")]|local,data = cax)